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ABSTRACT 

Upper limits on the shock speeds in supernova remnants can be combined with post-shock temperatures to 
obtain upper limits on the ratio of cosmic ray to gas pressure (Pcr/Pg) behind the shocks. We constrain shock 
speeds from proper motions and distance estimates, and we derive temperatures from X-ray spectra. The shock 
waves are observed as faint Ha filaments stretching around the Cygnus Loop supernova remnant in two epochs 
of the Palomar Observatory Sky Survey (POSS) separated by 39.1 years. We measured proper motions of 18 
non-radiative filaments and derived shock velocity limits based on a limit to the Cygnus Loop distance of 576 
± 61 pc given by Blair et al. for a background star. The PSPC instrument on-board /JOSAT observed the X-ray 
emission of the post-shock gas along the perimeter of the Cygnus Loop, and we measure post-shock electron 
temperature from spectral fits. Proper motions range from 2"7 to 5"4 over the POSS epochs and post-shock 
temperatures range from kT ^ 100-200 eV. Our analysis suggests a cosmic ray to post-shock gas pressure 
consistent with zero, and in some positions Pqr is formally smaller than zero. We conclude that the distance 
to the Cygnus Loop is close to the upper limit given by the distance to the backgroimd star and that either the 
electron temperatures are lower than those measured from ROSAT PSPC X-ray spectral fits or an additional 
heat input for the electrons, possibly due to thermal conduction, is required. 

Subject headings: ISM: individual (Cygnus Loop) - shock waves - supemova remnants - cosmic ray acceler- 
ation 



1. INTRODUCTION 

Supernova explosions are responsible for showering the 
surrounding interstellar medium (ISM) with heavy elements, 
influencing the distribution and composition of gas in the host 
galaxy. A typical supernova releases -^10^' ergs, sending a 
highly supersonic shock wave into its surroundings at speeds 
of up to 30,000 km s~' (Ghavamian et al. 2001). This ex- 
panding surge of energy heats and ionizes the neighboring 
gas out to hundreds of parsecs, initiating cloud collapse which 
leads to star formation and galactic evolution. Although a su- 
pernova is short-lived, the ejected remains propagate outward 
and interact with the ISM. These supernova remnants (SNRs) 
radiate over a large spectral range, transporting information 
about the blast wave morphology and local ISM. 

The Cygnus Loop is a well-known SNR whose distance is 
less than 576 ± 61 pc based on the distance to a star whose 
spectrum shows absorption features from shocked gas in the 
remnant (Blair et al., 2008). From a global perspective, its 
structure is governed by interactions with the surrounding in- 
homogeneous ISM. Despite its apparent spherical symmetry, 
both optical and X-ray observations suggest that the Cygnus 
Loop does not follow the simplified Sedov-Taylor morphol- 
ogy of an adiabatic blast wave expanding into a homogeneous, 
low density medium. Instead, the blast wave deceleration 
is the result of encounters with dense, extended clouds dis- 
tributed inhomogeneously in the pre-shocked gas (Levenson 
et al., 1997; Levenson et al., 1998). 

We focus on the non-radiative Ha shock filaments, primar- 
ily in the northeastern Cygnus Loop, several of which were 
previously observed by Hester et al. (1994), Blair etal. (1999, 
2008), Ghavamian et al. (2001) and Raymond et al. (2003). 

' Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138; salvesen@head.cfa.harvard.edu, jray- 
mond @ cfa.harvard.edu, edgar @ head.cfa.harvard.edu 

^ Department of Astronomy, The University of Michigan, 500 Church 
Street, Aim Arbor MI, 48109-1042; salvesen@umich.edu 



In a non-radiative shock, the cooling timescale of the hot 
shocked gas is much greater than the age of the shock. As the 
shocked gas has not had sufficient time to cool, post-shock 
radiation does not affect the shock evolution (McKee & Hol- 
lenbach, 1980). Therefore, non-radiative shock filaments pro- 
vide a laboratory for probing the conditions at the shock front 
before this information is lost by radiative or collisional pro- 
cesses. When a non-radiative shock encounters neutral gas, 
Balmer line emission is produced by the excitation of neutral 
hydrogen atoms prior to ionization. The resulting Ha fila- 
ments trace the outer edge of the Cygnus Loop non-radiative 
blast wave as it expands through low density, partiaUy neutral 
gas. These Ha filaments are useful for proper motion mea- 
surements of shock fronts and investigations of shocked gas 
parameters. 

SNR shock waves are believed to accelerate cosmic rays up 
to energies of 10^^ eV and to be responsible for a large frac- 
tion of interstellar cosmic rays, implying an efficient process 
for particle acceleration. The theory of diffusive shock accel- 
eration describes how charged particles (dominantly protons 
and electrons) passing back and forth across a shock front 
reach high energies, resulting in a power-law spectrum oc 
E"^ for cosmic rays above ~ lO' GeV (Blandford & Eichler, 
1987). This particle acceleration process requires a collision- 
less medium, as frequent collisions would return the velocity 
profile to Maxwellian. 

The efficiency of energy conversion in SNR shocks into 
the acceleration of cosmic rays is highly dependent on shock 
structure and is a topic of active research. Theory indicates 
that the fraction of the energy dissipated in a shock going 
into cosmic rays is bistable, being either very small or ^ 80% 
(Malkov et al. 2000). Recent calculations of the evolution of 
SNR shocks show that the efficiency starts off low and evolves 
graduaUy to somewhat smaller values than that predicted by 
steady shock models (Caprioli et al. 2008). Synchrotron 
emission in the radio and X-ray bands provides direct obser- 
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vation of particle acceleration in SNRs; however, the majority 
of the SNR energy going into particle acceleration is claimed 
by protons. Our current observational ability to detect 7-rays 
(and neutrinos) is limited, so we must settle for inferences of 
the total cosmic ray acceleration efficiency through indirect 
methods. Warren et al. (2005) inferred shock compression 
ratios greater than 4 from the distance between the blast wave 
and the contact discontinuity in Tycho's SNR. Their results 
suggest that the cosmic ray to gas pressure ratio, Pcr/Pg > 1) 
though projection effects might allow a smaller compression 
ratio (Cassam-Chenai et al. 2008). High cosmic ray accelera- 
tion efficiencies have also been inferred from the low electron 
temperature in 1E102-72.6 (Hughes et al. 2000) and from 
shock wave precursors to Bahner line filaments (Smith et al. 
1994; Hester et al. 1994). 

In this paper we find upper limits to Pcr/Pg by combin- 
ing measurements of post-shock electron temperature with 
the shock speed given by proper motions and distance upper 
limits. Our sample includes 18 non-radiative Ha shock seg- 
ments dominantly in the northeastern rim of the Cygnus Loop. 
This sample, although limited by the uncertainties in proper 
motion, temperature, and distance, provides constraints on 
Pcr/Pg- Several positions show apparently negative values 
of PcR, and we discuss the uncertainties in temperature and 
distance and the possibility of additional electron heating that 
could account for that non-physical result. 

2. OBSERVATIONS AND DATA REDUCTION 

We observe regions in the northern, eastern, and western 
Cygnus Loop over optical and X-ray bands. The Palomar Ob- 
servatory Sky Survery (POSS) obtained the optical data and 
the ROSAT Position Sensitive Proportional Counter (PSPC) 
obtained the X-ray data. We reduced two epochs of optical 
data for proper motion measurements. The X-ray data were 
obtained within two years of the more recent optical epoch. 
We matched WCS coordinates of all data to ensure consis- 
tency during optical and X-ray analysis. 

2.1. POSS Optical Data Reduction 

The POSS observed the Cygnus Loop over two epochs 
spanning 39.1 years. These photographic plates have been 
digitized by the Space Telescope Science Institute (STScI) 
and are publicly archived as the Digitized Sky Survey. We 
selected the filtered surveys whose transmission coefficient 
peaked near Ha (A = 6562.8 A), corresponding to surveys 
POSS-I E (red) and POSS-H F (red), to ensure detection of 
Ha shock filaments along the Cygnus Loop perimeter. We re- 
fer to the 1953 June 15 observations as POSS-1 and the 1992 
July 24 observations as POSS-11, with exposure times of 45.0 
and 120 minutes, respectively. The pixel scales of the dig- 
itized POSS-I and POSS-II plates are I'/l and I'/O, respec- 
tively. 

We analyzed proper motions of 18 Ha shock filaments 
which define the perimeter of the Cygnus Loop. We used IDL 
and SAOlmage DS9 for subsequent image analysis. From 
each POSS epoch, we hand selected a sub-image containing 
the filament of interest surrounded by reference background 
and foreground stars based on a common WCS (J2000) posi- 
tion. POSS-I and POSS-II pixel scales were resampled to 0"1 
using bilinear interpolation between grid points. Each set of 
sub-images was rotated appropriately to account for preces- 
sion effects during the time between observations. An addi- 
tional rotation was incorporated to aUgn the shock along the 
image axis for the proper motion analysis. Two-dimensional 



cross-correlation methods determined an optimal fine-tuned 
pixel shift to overlap the POSS-1 and POSS-11 sub-images. 
In the resulting images the background stars are stationary, 
leaving the propagating shock for analysis. All image manip- 
ulation was made relative to the POSS-I plate in an attempt to 
maintain the original integrity of the data. 

2.2. ROSAT PSPC X-Ray Data Reduction 

The Cygnus Loop X-ray data we analyzed were obtained 
with the PSPC detector on-board the /?0>SAr observatory. The 
observation timeline spans from 1992 Nov 04 to 1994 June 
03 with 3 separate pointings targeting 3 locations. Table 1 
lists the ROSAT PSPC observations relevant to our X-ray post- 
shock analysis of the Cygnus Loop. All observed data sets are 
publicly archived through the HEASARC. 

To prepare the ROSAT PSPC data for analysis, spectral re- 
gions were extracted from the PSPC events files requiring a 
minimum of 1,200 counts within each region. Source and 
background spectra were extracted and truncated to Pulse In- 
varient (PI) channels 0-255 using xselect version 2.4. We 
filter PSPC spectra to accept good events and, since we will 
be using statistics, we require a minimum of 20 counts per 
bin using the FTOOL grppha. We extracted X-ray emit- 
ting regions extending from 25" to 100" behind each Ha fil- 
ament of interest. As discussed below, we avoided the re- 
gion within 25" where the gas is significantly out of ionization 
equilibrium. Background spectra were taken outside the outer 
shock front and scaled by the ratio of the areas of the extrac- 
tion regions. PSPC response matrices were obtained from the 
legacy.gsfc.nasa.gov anonymous ftp archive and ancillary re- 
sponse files (ARFs) were generated using the FTOOL pcarf 
as outlined by Turner (1996). The resulting spectra were fit 
with XSPEC version 11.3.2 (Amaud, 1996). The PSPC de- 
tector is calibrated over the 0.1-2 keV soft X-ray energy band 
(Prieto, 1996). All spectra were fit in the 0.1-1.1 keV band 
as more energetic photons were not typically detected. The 
X-ray errors presented in this paper are at the 90% confidence 
level obtained with the XSPEC error command unless oth- 
erwise indicated. 

3. ANALYSIS AND RESULTS 

We surveyed 18 shock filaments, each of length 120", span- 
ning the eastern, northern, and western Cygnus Loop outer 
shell which is defined by its characteristic Ha emission. For 
each filament, we measured proper motion by comparing two 
epochs of POSS red digitized plates and post-shock gas tem- 
perature from spectral fits to ROSATPSPC X-ray data, assum- 
ing ionization equilibrium. Adopting the most recent distance 
upper limit to the Cygnus Loop of 576 ± 61 pc (Blair et al., 
2008), we constrain the ratio of cosmic ray to gas pressure. 
We discuss our methodology, the uncertainties in proper mo- 
tion, temperature and distance, and the results below. 

3.1. Proper Motion Measurements 

We determine the proper motion of 18 regions enclosing 
Ha shock filaments primarily along the northeastern perime- 
ter of the Cygnus Loop. Figures 1-6 show POSS-II images 
indicating the locations and dimensions of the extracted re- 
gions. In all images, north is at the top and east points to the 
left. Each extracted slit, outlined in blue with a corresponding 
identification number, has dimensions 25" x 120" and is po- 
sitioned parallel to a non-radiative Ha filament. Table 2 lists 
right ascension and declination relative to POSS-I for each of 
the 18 regions enclosing a non-radiative filament. Regions are 
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Fig. 1. — Image of the northern Cygnus Loop showing filaments 1-3. FiG. 4. — Image of the northeastern Cygnus Loop showing filament 14. 
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labeled by a filament number which increments with decreas- 
ing declination. We adopted the following Ha filament selec- 
tion criteria applicable to both POSS-I and POSS-II images: 

• Filaments must be bright relative to the local 
background. 

• Selected filaments must be continuous with no 
local branching structure. 

• Filaments must not precede another shock wave 
or strong optical emission within 100" to provide 
an extractable ROSAT PSPC post-shock region 
for X-ray temperature analysis. 

• A region enclosing a filament is not over- 
whelmed by stars. 

• Filament geometry appears to be rigid and pla- 
nar over both epochs allowing for uniform shock 
expansion in one dimension. 

Each filament is hand selected and both epochs of images 
are reduced as described in § 2.1 to align the shock along the 
image axis. Rectangular regions 25" x 120" parallel to and 
enclosing the filament in both epochs are selected for proper 
motion analysis. Within this region we extract slices which 
cross the shock but omit sub-regions containing stars. This 
method of avoiding stars ensures that the brightness of the 
propagating shock will dominate its surroundings in the se- 
lected sub-regions. Figure 7 shows a typical example of the 
region selection process just described. The black rectangu- 
lar regions of length 25" track shock propagation and omit 
stars. These black regions are stacked while regions enclosed 
by white bars are discarded to collect all starless sub-regions. 
A bootstrap method randomly selects half of the rows within 
these starless sub-regions to sample over the length of the fila- 
ment. For both POSS epochs, pixel values are summed down 
each column of sampled slices parallel to the shock, covering 
the 25" region width. A one-dimensional cross-correlation 
of the totaled rows between both POSS epochs determines an 
optimal offset corresponding to the shock propagation. The 
bootstrap process just described is repeated 1000 times to pro- 
vide a proper motion distribution and a statistical error esti- 
mate. The number of starless pixel rows did not exceed 600 
with a required minimum of 200, ensuring well over 1000 
possible combinations for bootstrap analysis. Figures 8 and 
9 show an example of a proper motion distribution after 1000 
trials and a plot of best-fit proper motions for a single iteration 
based on maximization of the correlation coefficient, respec- 
tively. 

Table 2 lists proper motion measurements for each filament 
with 90% confidence errors according to the standard devia- 
tion of the randomly sampled distribution (see Figure 8). The 
time elapsed between POSS-I and POSS-II observations is 
14,284 days (39.1 years). Proper motions over this timeline 
range from 2 "7 to 5 "4. We also give the conservative upper 
limit to the shock speed, v,, obtained from the upper limit to 
the proper motion (proper motion plus uncertainty) and the 
upper limit to the distance of the background star given by 
Blair et al. (2008) of 576H-61 = 637 pc. Matching filament 
locations (a, S) with corresponding proper motions in Fig- 
ures 1-6 demonstrates continuity of shock propagation along 
the rim of the Cygnus Loop and consistency of our measure- 
ments. Suitable filaments were restricted to the inner perime- 
ter of Ha emission to permit extraction of post-shock regions 
from ROSAT X-ray data. The post-shock gas must not 




Fig. 7. — Typical region selection process across an Ha filament. The 
shock is aligned along the image axis and enclosed by a 25" X 120" rectan- 
gular region. The black sub-regions are selected for analysis. The sub-regions 
enclosed by the white lines are dominated by stars; therefore, they are not ex- 
tracted for proper motion analysis. Identical sub-regions are extracted for 
both POSS epochs and the columns of the remaining starless regions are to- 
taled for cross-correlation comparisons. This selection process is specific to 
filament 7. 

encounter a trailing blast wave within ^ 100" to validate the 
assumption of ionization equilibrium required by our X-ray 
temperature analysis. While there are numerous, well-defined 
Ha filaments delineating the northeastern Cygnus loop, we 
only measure proper motions of those where the post-shock 
gas has adequate time to reach thermal and ionization equilib- 
rium, allowing us to constrain cosmic ray to gas pressure. 

3.1.1. Proper Motion Uncertainty 

The statistical uncertainties in the proper motions are de- 
termined by the measurement procedure to be 0." 1 to 0."2 for 
all but one of the filaments. Additional systematic uncertainty 
could arise from the cross registration of the images from the 
different epochs, though that error should be small because 
of the large number of stars used in the cross correlation. In- 
deed, we have compared the positions of individual stars in 
the different epochs and find them to differ by considerably 
less than an arcsecond. It is also possible that changes in 
brightness among different regions within an unresolved fil- 
amentary structure (see Blair et al. 1999) could cause errors 
in the proper motion. We expect such errors to be small, and 
measurement of 18 positions limits the impact of any such er- 
rors. Thus while there is a small systematic error, it is unim- 
portant because the uncertainties in temperature and distance 
dominate. 

We can also compare with other proper motion determi- 
nations. Shull and Hippelein (1991) reported proper mo- 
tions of 1"3 to 14."6 century"', with the latter value for a 
Balmer line filament close to our position 14, and Hester et 
al. (1986) found 7" per century for two Balmer line filaments 
in the northeast, one of which was observed by Blair et al. 
(1999). Blair et al. observed a non-radiative Ha filament in 
the Cygnus Loop (ayzooo = 20''56'«2'.7, 6j2om = 31°56'39"1) 
propagating into relatively dense neutral gas (see Figure 3). 
The filament proper motion was measured as 3 "6 ± 0"5 by 
comparing POSS-I and Hubble Space Telescope (HST) im- 
ages taken ^ 44.3 years apart using stars as stationary refer- 
ence points. However, the two reference stars used changed 
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Fig. 8. — Histogram distribution of proper motion measurements after 1000 
bootstrap iterations. Tlie bins are Of 1 wide and tlie results are normalized to 
unity. The 90% confidence range in proper motion is represented by a, given 
by multiplying the standard deviation by 1.645. This distribution is specific 
to filament 7. 
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Fig. 9. — Spread in correlation coefficient for one iteration between sam- 
pled POSS-I and POSS-n sub-regions contributing to the histogram in Figure 
8. Maximization of the correlation coefficient determines the optimal proper 
motion of the filament. This distribution is specific to filament 7. 

apparent position by 0"5 and O'.'l, increasing proper motion 
uncertainty. We measured this same filament with the method 
described in § 3.1 by combining POSS-I and two different 
POSS-II observations taken on 1992 July 24 and 1995 Au- 
gust 25 yielding proper motions of 2"6 ± 0"1 and 3"1 ± 
O'.'l, respectively. Asstiming a constant shock speed, we ex- 
trapolate these measurements to 3"0 and 3 "2 at the time of the 
1997 November 16 //5r observation. These two proper mo- 
tion measurements and that of Blair et al. (1999) were made 
relative to the POSS-I 1953 June 15 plate, providing three in- 
dependent measurements. Given that these results are within 
the limits presented in Blair et al. (1999) when reference star 
position changes are included, we consider our proper motion 
methodology justified. Although the //primage reveals intri- 
cate substructure, if we assume that the substructure of a par- 
ticular filament has not dramatically changed over the POSS 
epochs, the poor resolution of the POSS plates smoothes these 
features and should give accurate results. 



3.2. Post-Shock Temperature from ROSAT Spectral Fitting 

We extracted PSPC specfi-a of four 25" x 120" strips 
stacked behind and parallel to each filament for measure- 
ments of post-shock gas temperature behavior over 100". The 
single-temperature models apec and raymond show low 
post-shock temperature in the first ^ 25" behind each fila- 
ment. Within the 100" post-shock region, the fitted tempera- 
ture is consistentiy lower within the 25" region immediately 
behind the shock in agreement with Raymond et al. (2003). 
This suggests departure from ionization equilibrium immedi- 
ately following the shock because the gas has had insufficient 
time to ionize. 

We show below that the post-shock gas reaches ionization 
equilibrium after 25", so we extracted PSPC X-ray spectra 
of a 75" X 120" region extending from 25" to 100" behind 
each filament. Corresponding background regions were se- 
lected outside of the X-ray shock front and scaled, and source 
and background regions reduced as discussed above (§ 2.2). 
Figures 10-12 show the locations of extracted PSPC X-ray re- 
gions for all of the Hq filaments we considered, where north 
is at the top and east is toward the left. 

All of our spectral fits used the XSPEC model phabs to 
account for absorption of interstellar gas along the line of 
sight combined with a hot plasma X-ray emission model. The 
nominal interstellar column density value was Nh = 1.5 x 
10^" cm"^, but it was allowed to vary in all model fits (Ray- 
mond et al., 2003; Decourchelle et al., 1997). For all the soft 
X-ray fits presented in this paper, Nh is constrained within 

1 - 5 X lO^'^ cm"^. The kT^ and Nh parameters in the model 
are highly (and negatively) correlated, in the sense that a very 
small change in the temperature yields a large change in Nh- 
The fluctuations in A^^ are likely not significant. 

Ghavamian et al. (2000) measured Ha velocity line widths 
of a non-radiative shock located between filaments 5 and 6 
in the northeastern Cygnus Loop and found good electron-ion 
equilibration, T^/Ti ~ 0.7-1.0, in the post-shock region. We 
have assumed 7^ = 7] throughout our analysis. 

Single-temperature models and standard cosmic abundance 
assumptions are challenged by claims of bimodal temperature 
regions and depleted abundances (e.g. Nemes et al., 2008). 
Three fits were made to each PSPC X-ray spectrum: (i) a 
single-temperature fit assuming cosmic abundances (Anders 
& Grevesse, 1989), (ii) a single-temperature fit assuming de- 
pleted abundances, and (iii) a double-temperature fit assum- 
ing cosmic abundances. Depleted abundance fits fixed C, Mg, 
Al, Si, Ca, Fe, and Ni to 10% cosmic values on the basis 
that these elements are locked up in grains. Other abundances 
were fixed to cosmic values. 

Table 3 lists the resulting best-fit temperatures and 
values using the XSPEC models: 

(i) phabs X apec, 

(ii) phabs x vapec, 

(iii) phabs X (apec + apec), 

while Table 4 lists the same parameters using the XSPEC 
models: 

(i) phabs X raymond, 

(ii) phabs x vraymond, 

(iii) phabs X (raymond + raymond). 

Tables 3 and 4 compare the apec and raymond models 
to determine the sensitivity of the cosmic ray to gas pressure 
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ratio for slightly different X-ray temperature fits. Tcos and Tjep 
refer to temperature fits to models (i) and (ii), respectively. 
Tiow and Thigh are the bimodal temperatures from model (iii). 
Figures 13-15 show typical spectral fits for apec models (i), 
(ii), and (iii). 

Temperature estimates are required for calculations of the 
cosmic ray to gas pressure ratio downstream. The double- 
temperature models consistently produce better fits based on 
statistics, and Tiow remains continuous and consistent along 
ajacent filaments. However, there is no glaring physical rea- 
son to require two downstream temperatures. Hester et al. 
(1994) discuss how reflection of a blast wave from a dense 
cloud will enhance the temperature, and this may apply to the 
parts of the Cygnus Loop where the optical emission is bright, 
but we have avoided those complicated areas for this analysis. 
As discussed below, temperatures also depend on the assumed 
set of elemental abundances. We estimate the post-shock elec- 
tron temperature as r^os from the fits assuming cosmic abun- 
dances since they best represent the data. Lower limits on 
temperature are taken as Jiow from the double-temperature 
models, as this gives conservative lower limits to the gas 
pressure. Cosmic ray to gas pressure is inferred and an up- 
per limit is calculated based on these temperatures and mea- 
sured proper motions. Temperatures from both apec and 
raymond models are independently analyzed for all 18 post- 
shock regions as a consistency check on temperature measure- 
ments. 

3.2.1. X-Ray Temperature Uncertainty 

The derivation of temperatures from X-ray data, especially 
low-resolution X-ray data such as that provided by ROSAT, is 
fraught with systematic errors. In this section, we attempt to 
quantify some of those errors. 

Our extraction regions were chosen to provide at least 2,500 
good x-ray events in each, so that the statistical errors are min- 
imal. 

We have made two assumptions in our fitting: equilibrium 
ionization, and cosmic abundances. 

To investigate these, we produced spectra using the shock 
code of Cox and Raymond (1985) for a shock with = 
350 km s~\ which is typical of the velocities shown in Ta- 
ble 2, and a pre-shock density of 0.25 cm"^, which is typical 
of values determined from X-ray and UV observations (e.g., 
Raymond et al. 2003). We then integrated the surface bright- 
ness to distances which correspond to 25, 50, 75, and 100" 
behind the shock. The resulting spectra were multipled by a 
photoelectric absorption model with A^h = 1-5 x 10^" cm"^, 
convolved with the ROSAT response, and used with the 
XSPEC f akeit command to generate Poisson-sampled ran- 
dom spectra consistent with the shock model. We did this for 
both cosmic and depleted abundances. 

We then fit these spectra with a cosmic abundance, single 
temperature raymond model, assuming equilibrium ioniza- 
tion. We find that for the cosmic abundance cases, the fit- 
ted temperature underestimated the true input temperature by 
3%, which is a small effect. The first zone, within 25" of the 
shock, had a much lower fitted temperature. Accordingly, in 
our fitting of the Cygnus Loop data, we have ignored the first 
25" bin. 

Similarly, if we model a shock with depleted abundances, 
we find the spectrum is harder than with cosmic abundances. 
This is because much of the 1/4 keV band flux is provided by 
lines of Si, Mg, and Fe, while the 3/4 keV band flux at these 



shock speeds and temperatures is provided mostly by O and 
Ne which are undepleted even in dusty interstellar plasmas. 
Fitting such a spectrum with a model which presumes cos- 
mic abundances will overestimate the temperature to produce 
the harder spectrum. However, the model allows the column 
density of neutrals along the line of sight to vary, and this can 
compensate and even overcompensate for the excess emission 
in the soft band predicted by the cosmic abundance model, 
so that the temperature can be underestimated. Overall, we 
find that fitting a depleted shock spectrum to a cosmic abun- 
dance equilibrium model underestmates the temperature by 
about 11%. 

The large quantity of atomic physics data which go into 
simulations of the spectra of hot, collisional plasmas all have 
uncertainties as well. In order to assess these effects, we have 
fit the spectra with both the apec (Desai et al. 2005, Smith 
et al. 2001) and raymond (Raymond and Smith 1977, as up- 
dated) models. There is a systematic difference in the fitted 
temperatures of single-temperature, cosmic abundance mod- 
els, in the sense that the raymond model gives a fitted tem- 
perature about 0.15 ±0.05 keV lower than the apec model 
does. This is probably because apec includes only emis- 
sion lines for which reliable atomic data are available, while 
raymond includes emission lines which must be present, but 
for which the atomic data are more uncertain. This mainly 
affects the Mg, Si, S and Fe lines in the 1/4 keV band. 

Accordingly, we discuss lower limits to the electron tem- 
peratures. We use the low temperatures from the two- 
temperature, cosmic abundance fits as conservative lower lim- 
its to the temperature and therefore to the gas pressure. Even 
so in many cases we find that they imply upper limits to the 
cosmic ray pressure which are negative, a clearly unphysical 
situation. 

3.3. Other Temperature determinations for the Cygnus Loop 

The Cygnus Loop is one of the most studied supernova 
remnants in the Galaxy, and each X-ray observatory in turn 
has observed it. Each observatory has its own strengths and 
weaknesses, and the ensemble of results helps to constrain the 
systematics of any one observation. Many of these observa- 
tions cover the area we study, with special attention on the NE 
rim. 

Miyata et al. (2007) analyze Suzaku observations of a field 
in the NE. Their superior spectral resolution in the 0.3-2 keV 
band shows many lines of elements such as C, N, O, Ne and 
Mg. They fit two-temperature models to the spectra and find 
the higher temperature component, which is constrained by 
the lines in the spectrum, has kT greater than about 0.2 keV. 

Tsunemi et al. (2007) analyze data from a series of XMM- 
Newton pointings across the Cygnus Loop from NE to SW. 
They also fit two-temperature models to narrow concen- 
tric slices of the remnant. Their lower temperature (which 
they identify as the forward shock) is always greater than 
kT = 0.2 keV. 

The Chandra observatory has also observed the NE rim. 
Katsuda et al. (2008) report fits to the XSPEC vpshock 
model, which allows for nonequilibrium ionization effects, 
and variable abundances. They also obtain fitted temperatures 
uniformly in excess of 0.2 keV. 

It therefore seems that these high temperatures are quite 
robust features of the X-ray emission from the NE rim of 
the Cygnus Loop. Chandra's high spatial resolution, XMM- 
Newton's large collecting area and Suzaku's superior spectral 
resolution aU lead to the same conclusion we obtain from the 



Cygnus Loop 



7 




'■ 1 


^ 

^^^^ - 




y 1 1 1 1 \ — 1 — h- 

1 ■ 



Fig. 10. — ROSATPSPC image of the northeastern Cygnus Loop showing 
locations of the X-ray source extraction regions behind filaments 1-14. 



0.2 0.5 1 

Energy (keV) 

Fig. 13. — Representative example of a ROSATPSPC single-temperature 
apec model fit with abundances fixed to cosmic values. In most cases, this 
model was statistically better than the depleted model. This fit is specific to 
filament 7. 
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Fig. 11. — ROSAT PSPC image of the western Cygnus Loop showing lo- 
cations of the X-ray source extraction regions behind filaments 15 and 16. 



Energy (keV) 

Fig. 14. — Representative example of a ROSAT PSPC single-temperature 
apec model fit with C, Mg, Al, Si, Ca, Fe, and Ni abundances depleted to 
10% cosmic values and fixed. The rational for depletion is that these elements 
may be locked up in grains. This fit is specific to filament 7. 




Fig. 12. — ROSAT PSPC image of the eastern Cygnus Loop showing loca- 
tions of the X-ray source extraction regions behind filaments 17 and 18. 
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Fig .15. — Representative example of a ROSAT PSPC double-temperature 
apec model fit with abundances fixed to cosmic values. In most cases, a 
bimodal temperature distribution was statistically significant, providing mo- 
tivation to investigate this in future works. This fit is specific to filament 7. 
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ROSAT data: kT values of 150 eV or higher are ubiquitous in 
the non-radiative filaments of the Cygnus Loop. 

3.4. Distance 

The distance to the Cygnus Loop is a substantial uncer- 
tainty, and our measurements can be used to obtain a lower 
Umit to the distance. Minkowski (1958) obtained the canoni- 
cal value of 770 pc from the velocity ellipse and proper mo- 
tions, though Braun and Strom (1986) used the same data to 
obtain 460±160 pc using the mean expansion velocity instead 
of the extreme. Shull & Hippelein (1991) obtained a distance 
of 600 pc based on Fabry-Perot scans of a large number of 
fields in the Cygnus Loop, but with a range of 300 to 1200 pc. 
Sakhibov and Smimov (1983) also compared proper motions 
and radial velocities to obtain an estimate of 1400 pc. Though 
the method of comparing radial velocities and proper motions 
relies only upon the assumptions of symmetric expansion and 
radial motion, it leads to a wide range of distance estimates. 

An independent distance estimate by Blair et al. (1999) 
combined the proper motion of a single filament with a shock 
speed obtained from spectroscopic analysis and fits to shock 
models. They derived a distance of 440 pc with a range of 340 
to 570 pc. Their method assumed that Pcr/Pg is zero, and a 
larger value would imply a higher shock speed and therefore a 
larger distance. This proper motion had a relatively large un- 
certainty because of the possible motions of the two reference 
stars. 

The most solid limit on the distance to the Cygnus Loop 
comes from Blair et al. (2008). They used FUSE to observe 
an sdO star and found strong high velocity O VI absorption 
lines matching the O VI emission in the adjacent nebulosity. 
A fit to the FUSE spectrum provided the temperature and sur- 
face gravity of the star, indicating a distance of 576±61 pc. 
The upper end of this range, 637 pc, should be a firm up- 
per limit to the distance of the Cygnus Loop. Distances to 
sdO stars are generally difficult to determine, but the spectral 
fits should give an accurate temperature. The uncertainty in 
the surface gravity probably dominates the 10% uncertainty 
in the distance. 

3.5. Analyses of other Balmer line shocks in the Cygnus Loop 

Several of the Balmer line shocks at the periphery of the 
Cygnus Loop have been studied previously. The shock ob- 
served by Blair et al. has been observed extensively at optical 
and UV wavelengths (e.g. Raymond et al. 1983; Fesen & 
Itoh 1985; Long et al. 1992; Hester et al. 1994). We do not 
include it here (except to check our proper motion measure- 
ments; § 3.1.1) because it is too slow to produce X-ray emis- 
sion. However, it fits in with our analysis in that the shock 
speed is estimated to be 150 to 190 km s~' and the proper mo- 
tion is around two thirds the typical values we measure for the 
X-ray producing shocks. 

The filament corresponding to our positions 15 and 16 was 
studied by Raymond et al. (1980) and Treffers (1981), who 
measured an Ha line width corresponding to a shock speed 
faster than 170 km s"' . The filament corresponding to our po- 
sitions 17 and 18 was discussed by Fesen, Kwitter & Downes 
(1992) and Graham et al. (1995). The latter paper finds that 
the Balmer line filament in the region we observed is indeed 
the blast wave, rather than part of the shock structure refracted 
around the dense cloud just to the north. 

Detailed studies of a position in the northeast have been 
made by Ghavamian et al. (2001) and Raymond et al. (2003). 



We have not included this exact region because two filaments 
overlap near there, but it lies between our regions 5 and 6. 
Ghavamian et al. obtained Ha and H/? profiles and found 
nearly equal proton and electron temperatures {Te/Ti > 0.7 
and a shock speed of 235 to 395 km s"'. Raymond et al. 
(2003) combined those results with UV spectra from FUSE 
and the Hopkins Ultraviolet Telescope. They found further 
evidence for near equilibration among the particle species, 
with an oxygen kinetic temperature near that of hydrogen 
ijo < 2.5Tp). They also analyzed the ROSAT X-ray spectrum, 
finding 0.14 < Te <0.2 keV. They also observed a bright rim 
about 25"thick which they attributed to enhanced emission 
in the 0.25 keV band due to non-equilibrium ionization. They 
found that a 350 km s"' shock matches the optical, UV and X- 
ray data. This speed combined with the average of the proper 
motions of positions 5 and 6, 3f'75, would give a distance of 
750 pc. It should be mentioned, however, that van Adelsberg 
et al. (2008) computed models of Ha profiles in non-radiative 
shocks and were unable to simultaneously match the broad 
Une component width and the broad to narrow intensity ratio 
reported by Ghavamian et al. 

4. COSMIC RAY TO GAS PRESSURE 

The proper motion and temperature measurements de- 
scribed above, combined with a distance to the Cygnus Loop, 
allow calculation of the upper limit to the cosmic ray to gas 
pressure ratio downstream. 

Taking the perspective of the shock front, momentum con- 
servation applies to the pre- and post-shocked gas: 



d\ 

p— +p(v • V)v+ VP+ pVt/i = 
at 

Assuming a steady solution exists, ignoring the gravita- 
tional force term resulting from a field, 0, and applying the 
continuity equation pv = constant, the momentiun equation 
becomes, 

dv dP d , 9, „ 
pv— + —=—(P+pv^) = Q. 
dx dx dx 

This translates into a statement of constant momentum flux 
across the shock, allowing us to equate pressure terms, 

P\ + ^ + P\v] = Pcr + Pg+^ + P2vI, 
87r 87r 

where Vj is shock speed, p is density, Pqr is cosmic ray pres- 
sure originating in a dissipative shock envirormient, Pq is 
downstream gas pressure, B^/Sw is magnetic pressure, and the 
subscripts 1 and 2 represent conditions in the pre- and post- 
shock regions, respectively. We assume the ram and down- 
stream pressures dominate, neglecting the ambient gas pres- 
sure and interstellar magnetic field. Assuming the case of a 
strong shock with adiabatic exponent 7 = ^ where 112/111 =4 
and V2 = \vs, along with the equations of state pi = nim and 
Pg = «2^72. we derive: 



l-h 



3mVj 
l6kT2 



'imp^d^ 



where T2 is the mean post-shock temperature. We assume ion- 
ization equiUbrium exists, and that = T^. = Ti (Ghavamian 
et al., 2000). We calculate shock velocities assuming oiily 
a tangential component, = = where ^ is the angu- 
lar propagation across the sky over a particular timescale (i.e. 
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proper motion) and d is distance to the Cygnus Loop. Adopt- 
ing cosmic abundances with n(He)/n(H) = 0.098 (Anders & 
Grevesse, 1989), the mean mass per particle behind the shock, 
m, is determined by dividing the weighted molar mass per nu- 
cleus by the average number of particles per nucleus. The 
cosmic ray to gas pressure ratios presented in this paper are 
computed according to the equation derived above. We note 
that the compression factor of 4 may be an underestimate if 
cosmic ray pressure significantly impacts the gas dynamics, 
decreasing the post-shock flow speed relative to the shock, 
but as will be seen, we are primarily interested in low effi- 
ciency shocks. We neglect line of sight velocity components 
in our calculations, assuming the shock velocity is mostly 
perpendicular to the line of sight. We note that the appear- 
ance of multiple propagating shock layers is an artifact of our 
two-dimensional perspective. Instead, the complicated three- 
dimensional structure of the Cygnus Loop is wavy, warped, 
and sheet-like (Raymond et al. 2003). Cosmic ray to gas pres- 
sure ratios were calculated by combining proper motion and 
temperature measurements with an adopted distance of 576 ± 
61 pc, obtained from a distance upper limit to a background 
subdwarf OB star (Blair et al., 2008). 

Note that the density has dropped out of the equations. The 
results we present here do not depend upon the ambient den- 
sity, which is not very well known. 

Table 2 lists best-fit values and upper limits of Pcr/Pc for 
all 18 filaments calculated independently with the apec and 
raymond ionization equilibrium temperature fits. To obtain 
conservative upper limits, we have used the upper limits to 
the shock speed v, from Table 2 based on upper limits to 
the proper motions and the upper limit to the distance of 637 
pc (Blair et al. 2008). We find several negative values for 
Pcr/Pc- This result is unphysical and requires further inves- 
tigation. For all instances of negative pressure ratios we set 
Pcr/Pc to zero and derive a minimum distance estimate to 
the Cygnus Loop based on our proper motion and tempera- 
ture uncertainties. 

5. DISCUSSION 

Inspection of the derived cosmic ray to gas pressure equa- 
tion emphasizes the importance of constraining proper mo- 
tion, distance, and temperature to infer tight upper limits on 
Pcr/Pc- Pcr/Pc scales with the square of the shock speed, 
given by the product of proper motion and distance. In gen- 
eral, the upper Umits indicate a small value for the ratio, and 
in some cases it is formally negative. 

In the context of the discussions of the uncertainties given 
above, it is clear that the proper motion measurements cannot 
account for the unphysical values of Pcr/Pc- Most distance 
estimates are significantly smaller than the 637 pc we have 
adopted as the nominal upper Hmit, and it seems like a strain 
to make the distance large enough (^^ 1 kpc) to make aU the 
upper limits to Pcr/Pc positive. However, it may be possible 
given the general uncertainties in SdO star distances. 

Therefore, we conclude that the heart of uncertainty likely 
lies in post-shock temperature measurements and electron-ion 
thermal equilibration assumptions. Assuming ionization equi- 
librium and that shock velocity measurements are reasonably 
accurate, we require lower post-shock electron temperatures 
than measured from out ROSATPSPC X-ray spectral fits. The 
fits based on the Raymond & Smith (1977) code give temper- 
atures about 10% lower than those using APEC, but that is 
several times too smaU to explain the discepancy. 



5.1. Heating by Thermal Conduction 

There is another possible way out of the discrepancy. We 

have assumed that Tg = Tj behind the shock, and it is possible 
that Te > Ti. However, such behavior is not seen in shocks in 
the interplanetary medium, and there is no obvious physical 
reason for such electron heating in the shock. Ghavamian et 
al. (2006) describe a wave heating mechanism that would 
produce particle temperatures somewhat higher than observed 
here, but Tp from Ghavamian et al. (2001) and a higher 
would then require a large value of Vs. 

An interesting alternative is electron heating by thermal 
conduction from gas farther behind the shock. A number of 
papers have explored the global structure of SNRs with ther- 
mal conduction (e.g., Slavin & Cox 1992; Cox et al. 1999; 
Shelton et al. 1999). We can obtain an observational esti- 
mate from X-ray observations. Nemes et al. (2008) show 
the variation of Te with radius in the northeastem Cygnus 
Loop based on XMM spectra, with a gradient of 0.05 keV 
over a distance of 0.07 times the shock radius, or WT ~ 
1.8 X 10"'^ K cm"'. This implies a volumetric heating rate 
of about 3.4 X 10"^^ erg cm"^ s"', which could heat the gas 
by about 3 x 10^ K as it travels over the 100" region we an- 
alyze. Thus the true post-shock temperature could be some- 
what lower than the X-ray temperatures we measure by ap- 
proximately the amount required for Pcr greater than zero. 
The above estimate requires a temperature gradient parallel 
to the magnetic field, so it could not be correct everywhere. 
However, only 7 of the 18 positions show Pcr nominally 
smaUer than zero. 

6. SUMMARY 

Constraining Pcr/Pc in supernova shocks is important for 
assessing the efficiency of energy dissipated by the SNR into 
accelerating cosmic rays. We combine measured proper mo- 
tions with temperatures derived from X-ray spectra and an up- 
per limit to the distance to obtain upper limits to Pcr/Pc- We 
measured proper motions and post-shock temperatures of 18 
faint nonradiative Ha filaments in the Cygnus Loop. Proper 
motion measurements of filaments based on image matching 
and correlation techniques yield continuous results over ex- 
tended shock regions. Post-shock electron temperatures from 
X-ray fits are higher than expected if the current 576 ± 61 pc 
distance measurement to the Cygnus Loop is accurate. 

We conclude that 1) Pcr/Pc is small, 2) the distance to the 
Cygnus Loop must be close to or even larger than the upper 
limit given by the apparent distance to a star that lies behind 
the SNR, 3) uncertainties in the temperature derived from X- 
ray spectra might dominate the uncertainties in the analysis, 
and 4) thermal conduction from hotter interior gas could alter 
the immediate post-shock temperature enough to account for 
the observations. 

Future work will focus on constraining post-shock temper- 
ature and investigating the validity of ionization equilibrium 
on a more global scale along the northeastern Cygnus Loop, 
and placing tighter constraints on the distance to the Cygnus 
Loop. 
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Fig. 16. — Mapping filaments by their numbers defined in Figures 1-6 corresponding to Pcr/Pg best fits with error bars extending to the calculated upper limit. 
Results for XSPEC model fits apec are shown as filled symbols and raymond results are shown as open symbols. Groupings of local filament regions cluster 
around similar Pcr/Pg values. Filaments 4-13, which delineate the northeastern rim of the Cygnus Loop, appear to show Pcr/Pg values consistent with zero. 



TABLE 1 



Cygnus Loop ROSAT PSPC Observation Parameters 



Name 


«y20oo 


5y20oo 


Date 


Exposure 








(UT) 


(ksec) 


rp500034a01 


20 55 02.4 


-H32 07 12.0 


1992 Nov 04 


12.0 


rp500267n00 


20 56 21.6 


-H30 24 00.0 


1993 Nov 14 


9.24 


rp500268n00 


20 45 40.8 


-H31 02 24.0 


1994 Jun 03 


9.75 



Note. — The coordinates listed describe the target location for a particular i?05Ar PSPC pointing. Units of right ascension are hours, minutes, and seconds. 
Units of declination are degrees, arcminutes, and arcseconds. We adopt this system of units ttu-oughout this paper. 
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TABLE 2 



Proper Motion and Pcr/Pg Values for Selected Ha Filaments 



Filament ID 




5y2000 


Proper Motion 
arcsec / 39.1 yr 


km s ' 


(PcR/PG)m>ec 


(Pcr/Pg)S? 


jMIN 
apec 

pc 


(PrR / /V;)rflvninnH 


(Pcr/Pg)"'^^ a 


jMIN 
raymond 

pc 


1 


20 51 23.9 


32 24 22.5 


S'.'l ±0f'2 


403 


-0.032 


0.408 


537 


0.111 


1.274 


422 


2 


20 51 29.9 


32 24 21.8 


5f'4 ± Q'.'l 


433 


0.091 


0.625 


500 


0.190 


0.637 


498 


3 


20 51 38.6 


32 24 14.5 


5."2 ± 0."2 


416 


-0.032 


0.501 


520 


0.053 


0.582 


507 


4 


20 54 13.1 


32 21 14.0 


2"7 ± 0."2 


225 


-0.723 


-0.616 


1028 


-0.698 


-0.560 


960 


5 


20 54 26.2 


32 19 36.5 


3^'4 ± 0^'2 


278 


-0.543 


-0.344 


786 


-0.473 


-0.325 


775 


6 


20 54 43.4 


32 16 04.0 


A'.'l ± Q'.'l 


333 


-0.337 


-0.004 


638 


-0.239 


0.550 


512 


7 


20 55 06.3 


32 10 03.8 


3'.' a ± Q'.'l 


240 


-0.587 


-0.496 


897 


-0.561 


-0.193 


709 


8 


20 55 14.7 


32 07 38.9 


3f'4±0f'l 


274 


-0.507 


-0.428 


843 


-0.456 


-0.397 


821 


9 


20 55 18.9 


32 06 58.0 


3'.'1 ± Q'.'l 


294 


-0.384 


-0.275 


748 


-0.338 


-0.259 


740 


10 


20 55 34.6 


32 01 40.2 


3'.'5±Q'.'l 


279 


-0.400 


-0.317 


771 


-0.383 


0.026 


629 


11 


20 55 44.9 


31 59 43.8 


Z'.'l ±Q'.'l 


254 


-0.507 


-0.432 


845 


-0.424 


0.066 


617 


12 


20 55 51.8 


31 57 34.8 


A'.'Q ± Q'.'l 


319 


-0.231 


-0.115 


677 


-0.203 


0.253 


569 


13 


20 55 56.5 


31 55 55.1 


4:'2 ± o:'7 


375 


-0.120 


0.256 


568 


0.020 


1.034 


447 


14 


20 57 20.2 


31 37 32.4 


4'.' 3 ± 0."1 


342 


-0.176 


0.030 


628 


-0.127 


0.299 


559 


15 


20 45 11.9 


31 03 50.0 


3^'4±0^'l 


272 


-0.415 


-0.351 


791 


-0.329 


-0.182 


704 


16 


2045 15.3 


31 01 41.4 


3'.'3±Q'.'l 


264 


-0.492 


-0.443 


854 


-0.470 


-0.428 


842 


17 


20 56 37.4 


30 08 34.4 


4f'5 itOf'l 


358 


-0.054 


0.049 


622 


0.013 


0.434 


532 


18 


20 56 34.8 


30 06 27.8 


4-'.' &± Q'.'l 


379 


-0.044 


0.237 


573 


0.099 


0.707 





Note. — The coordinates listed represent the right ascension and declination at the center of the extracted filament. The measured proper motions are derived 
from comparing POSS-1 and POSS-II images observed 39.1 years apart. Errors on proper motion do not include a < Of'l uncertainty from image alignment. 
Shock speed, v,, is calculated from the product of proper motion and distance using the upper limits of proper motion + uncertainty and 576+61 = 637 pc (Blair 
et al. 2008). All cosmic ray to gas pressure ratio calculations are based on these conservative upper limits. Minimum distances to the Cygnus Loop based on our 
measurements assume Pcr/Pg = with an upper limit on proper motion and lower limit on temperature. We compare results from individual X-ray temperature 
fits with the XSPEC models apec and raymond. 
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TABLE 3 



X-ray Spectral Fit Parameters with the ape c Model 



Filament ID 


T 

A cos 

eV 


Mil 

lO^" cm-2 


VA. / ^ /COS 


^oep 

eV 




eV 


Thigh 

eV 


Va / '^ydoubie 


1 


187!« 


2 440.2 


1.65 


162^10 


3.37 


137 


536 


1.13 


2 


192^7 


2 4+0.2 
^•^-0.2 


1.80 


198!io 


3.01 


137 


602 


1.28 


3 


198t^ 


2540.2 
^•-'-0.2 


1.44 


230+;o 


2.47 


137 


336 


1.16 


4 


190!« 


21!40.2 
^•-'-0.2 


0.96 


194!| 


1.86 


157 


454 


0.85 


5 


18l!i 


2 5+<'-2 
^■-'-0.2 


1.48 


159!^ 


3.21 


140 


605 


1.06 


6 


179!^ 




2.33 


151 


5.75 


132 


380 


1.60 


7 


154!| 




2.22 


119 


5.95 


136 


541 


1.53 


8 




1.2 


1.26 


118!^ 


3.44 


156 


661 


1.06 


9 


159^f 


^•'-0.2 


1.45 


115!^ 


3.63 


141 


656 


1.14 


10 


146 


2.1 


5.59 


80.8 


10.7 


135 




3.90 


11 


139 


1.9 


5.07 


80.8 


9.92 


135 


690 


3.95 


12 


145« 




2.91 


80.8 


7.18 


136 


685 


1.55 


13 


139 


1.7 


5.03 


80.8 


9.34 


133 


863 


3.40 


14 


159+8 


1 8+<'-3 

'■°-0.3 


1.48 




3.14 


135 


727 


0.89 


15 


140^3 


2 1+0-2 
^■'--0.2 


2.81 


80.8 


5.24 


135 


812 


2.29 


16 


150:^ 


2 4+0.1 


1.02 


86.8 


1.91 


148 


692 


1.03 


17 


152!" 




1.19 


146!« 


1.91 


145 


159 


1.24 


18 






1.32 




1.88 


138 


633 


1.24 



Note. — ROSAT PSPC X-ray spectral parameters corresponding to fits behind each filament with the XSPEC model apec. The subscript 'cos' refers to 
single-temperature model fits phabs X apec with abundances fixed to cosmic. The subscript 'dep' refers to single-temperature model fits phabs X vapec 
with depleted abundances fixed to 10% cosmic. The subscript 'double' refers to double-temperature model fits phabs X (apec+apec) with abundances 
fixed to cosmic. We allowed Nh to vary from an initial value of 1.5 X 10^" cm"^ Wg stress that the errors listed in the table are generated in XSPEC and are 
not representative of actual uncertainties based on the variation in best fit values when comparing similar models. Fits with no errors or parameters listed are 
unphysical or limited by statistics. 
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TABLE 4 



X-ray Spectral Fit Parameters with the raymond Model 



nentID 


Tcos 

C V 


Nh.cos 
jq20 




Tdep 

C V 


(X^/f')dep 


C V 


Thigh 

C V 


(x7^)dc 


1 
1 




T c+0.5 
-'•-'-0.4 


i.oy 


/I c -3+22 
^-'•^-20 


1 ^7 








z 


17(^+8 

—3 


-1 n+0.2 




405+27 


1. /U 


1 jD 


1(Y) 
/VIZ 


1 97 
l.Z / 


J 


—5 


Q+0.2 


1 Al 


407+26 


z. 01 






1 1 ^ 


4 


174!«j 


'I o+fl fi 

2.8!^:° 


1.03 


AAA +20 

AAA T^-^v 

**^-16 


1.93 


137 


448 


0.84 


J 


i J t _^ 


-1 fi+0.4 


1. J i 




9 IQ 

Z. I7 


t "Xfii 
ljo 




1 nft 


D 




-1 9+0.2 
•'■^-0.2 


9 1 7 


477+9 


1 1 8 


84 8 


266 


1 28 


7 


it J 2 


T c+0.2 


9 (19 
z.uz 


079+12 
J//_(2 


2 85 


84 9 


292 


1 no 


Q 
O 




9 1+0.2 


1. iO 




9 90 

z.zy 


1 




1 ns 

l.Uo 


Q 




T c+0.2 
^•-'-0.2 


1. J'f 


41 S+17 
413_i7 


9 09 

z.uz 


I JO 




111 

1. i 1 


1 n 


1 AO 


z.o 


S 10 

J. ly 


■J7-3+10 


9^^ 


SO 


jjy 


9 ftO 
Z.DU 


1 1 
i 1 


1 10+5 

iiy_5 


q £+0.4 




"^^fi+S 


9 '^9 
Z. JZ 


71 


AAQ 


1 ft7 
l.D / 


12 


140^3 


1 Q+0.1 

^•^-0.2 


2.58 


379^10 


2.08 


96.1 


425 


1.04 


13 


120+5 


11 ■1+0.4 
^•"^-0.3 


4.63 


366:li 


2.98 


82.1 


627 


1.84 


14 


150!f 


^•°-0.3 


1.34 




1.17 


107 


670 


0.84 


15 


122!« 


o 7+0.6 
•^■'-0.5 


2.66 


114!^ 


2.83 


107 


704 


1.81 


16 


144+^ 




1.02 


133!| 


1.16 


144 


145 


1.06 


17 


142:11 


5 3-^-5 

-^•^-0.6 


1.22 


143^ 


1.39 


106 


693 


1.03 


18 


147!io 


-^•^-0.8 


1.32 


150ti§ 


1.52 




147 


1.39 



Note. — ROSAT PSPC X-ray spectral parameters corresponding to fits behind each filament with the XSPEC model raymond. The subscript 'cos' refers 
to single-temperature model fits phabs X raymond with abundances fixed to cosmic. The subscript 'dep' refers to single-temperature model fits phabs X 
vraymond with depleted abundances fixed to 10% cosmic. The subscript 'double' refers to double-temperature model fits phabs X (raymond+raymond) 
with abundances fixed to cosmic. We allowed Na to vary from an initial value of 1.5 X 10^" cm"^. We stress that the errors listed in the table are generated 
in XSPEC and are not representative of actual uncertainties based on the variation in best fit values when comparing similar models. Fits with no errors or 
parameters listed are unphysical or limited by statistics. 



